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We present a two-dimensional classical stochastic differential equation for a displacement field of 
a point particle in two dimensions and show that its components define real and imaginary parts of a 
complex field satisfying the Schrodinger equation of a harmonic oscillator. In this way we derive the 
discrete oscillator spectrum from classical dynamics. The model is then generalized to an arbitrary 
potential. This opens up the possibility of efficiently simulating quantum computers with the help 
of classical systems. 
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1. In a recent analysis of quantum mechanics from the 
point of view of information processing, one of us ||] 
pointed out that decoherence will become an insurmount- 
able obstacle for the practical construction of quantum 
computers. It was suggested that instead of relying 
on quantum behavior of microparticles it seems more 
promising to simulate quantum behavior with the help 
of fast classical systems. As a step towards such a goal 
we construct a simple classical model which allows us to 
simulate the quantum behavior of a harmonic oscillator. 
In particular we show that the discrete energy spectrum 
with a definite ground state energy can be obtained in 
a classical model. In the latter respect we go beyond an 
earlier model in Ref. |^ whose spectrum had the defect 
of being unbounded from below. Finally, the model is 
generalized to an arbitrary potential. 

2. For a point particle in two dimensions we de- 
fine a time-independent displacement field u(x) — 
(it^(x), u^(x)) parametrized by the spatial coordinates 
X — {x^,x'^). The reparametrization freedom is fixed by 
choosing harmonic coordinates in which 



Au(x) = 0, 



(1) 



where A is the Laplace operator. This condition implies 
that the components m^(x) and u^(x) satisfy the Cauchy- 
Riemann equations 



{fi,i^,... = l,2), 



(2) 



where e^i/ is the antisymmetric Levi-Civita pseudoten- 
sor. The metric is 5^^, so that indices can be sub- or 
superscripts. 

The particle is supposed to be in contact with a heat 
bath of "temperature" h. Its classical orbits x(t) are 
assumed to follow a stochastic differential equation con- 
sisting of a fixed rotation and a random translation in 
the diagonal direction n = (1,1): 



where iv is the rotation vector of length uj pointing or- 
thogonal to the plane, and r](t) a white-noise variable 
with zero expectation and the correlation function 



(4) 



For a particle orbit x(t) starting at x(0) = x, the position 
x(t) at a later time i is a function of x and a functional 
of the the noise variable ri{t') for < t' < t: 



x(t) = Xt[x;?7] 



(5) 



To simplify the notation we indicate the time dependence 
of functionals of 77 by a subscript t. 

We now introduce a time-dependent displacement field 
u(x; t) which at t = is equal to u(x) and evolves with 
time as follows: 

u(x;i) = ut[x;77] = u(Xt[x;77]) , (6) 

where the notation [x; rj] indicates the variables as in 
©• 

As a consequence of the dynamic equation (^), the 
change of the displacement field in a small time inter- 
val from t = to t — At has the expansion 



Auo[x;77] ^ At{[LO xx]- V}uo[x;f]] 
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dt'rj{t') (n- V)uo[x;77] 
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The omitted terms are of order At'^/^ . 

We now perform the noise average of Eq. (Q) , defining 
the average displacement field 
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u(x;<) = (ut[x;r7]). 



(8) 



AS(x) = 0, 



(18) 



Using the vanishing average of ri{t) and the correlation 
function (Q), we obtain in the hmit A< ^ the time 
derivative 



5tu(x;0) = Hu(x;0) 
with the time displacement operator 



H^{[a;xx].V} + -(n.V)2 



(9) 
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The average displacement field u(x; t) at an arbitrary 
time t is obtained by the operation 



u(x; t) = U{t)u{x; 0) = e^*u(x; 0) 



(11) 



Note that the average over rj makes the operator 7i time- 
independent: li.U{t) = b{{t)'H. Moreover, the operator 
H commutes with the Laplace operator A, thus ensuring 
that the harmonic property ([^) of u(x) remains true for 
all times, i.e., 



Au(x;t) = 0, 



(12) 



3. We now show that Eq. (||) describes the quantum 
mechanics of a harmonic oscillator Let us restrict our at- 
tention to the line with arbitrary = x and X2 — 0. Ap- 
plying the Cauchy-Riemann equations (^, we can rewrite 
Eq. (w in the pure a;-form 



(13) 
(14) 



Now we introduce a complex field 

where we have written uj' (x) for (x) \xi=x,x2=o- This 
satisfies the differential equation 



ihdtip{x; t) 



^dl + —x^-—]i^{x;t), (16) 



which is the Schrodinger equation of a harmonic oscilla- 
tor with the discrete energy spectrum En = {n+ l/2)huj, 
n = 0,1,2,... . 

5. The method can easily be generalized to an arbitrary 
potential. We simply replace (||) by 



xHt) = ~d2S\x{t))+n^vit), 



(17) 



where S(x) shares with u(x) the harmonic property (j^): 



i.e., the functions 5^(x) with /i = 1,2 fulfill Cauchy- 
Riemann equations like m''(x) in (||). Repeating the 
above steps we find, instead of the operator (fol). 



n = ~{d2S')d, - id,s')d2 + -(n • v)2, 

and Eqs. ( [T^ ) and ( |l4|) become: 



dtul^ {dxS^)dxu1 ~ -dlu^ 



(19) 

(20) 
(21) 



This time evolution preserves the harmonic nature of 
u(x). Indeed, using the harmonic property AS(x) = 
we can easily derive the following time dependence of the 
Cauchy-Riemann combinations in Eq. (^): 

dt{diui - d2U2) = H{diui - 92^2) 

- d2diS\diUi - d2U2) + diS\d2Ui + diU2), 

dt{d2Ui + diU2) = ii{d2Ui + diU2) 

- d2dlS\d2Ui + dlU2) - diS\dlUi - d2U2). (22) 

Thus diUi — d2U2 and d2Ui + diU2 which are zero at any 
time remain zero for all times. 

On account of Eqs. (pT]), the combination 



il)(x;t) 



[U^ {x;t) + iu'^ {x;t)] 



satisfies the Schrodinger equation 



ifidtip{x; t) = 



-^dl + V{x) 



ip{x;t), 



(23) 



(24) 



where the potential is related to S^{x) by the Riccati 
differential equation 



V{x)^l[dxS\x)r-^dlSHx). 



(25) 



The harmonic oscillator is recovered for the pair of 
fimctions S'i(x) + iS^{x) = lu{x^ + ix^ f 12. 

5. The noise r^it) in the stochastic differential equation 
Eq. ( |l7| ) can also be replaced by a source composed of 
deterministic classical oscillators g/c(i), fc = 1, 2, . . . with 
the equations of motion 



Qk^Pk, Pk 



as 



(26) 



(27) 
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The initial positions qk{0) and momenta Pfc(O) are as- 
sumed to be randomly distributed with a Boltzmann fac- 
tor e-l3Husc/n^ g^pj^ ^Yia,t 

(gfe(O)gfe(O)) = h/iul (pfe(OK(0)) = h. (28) 

Using the equation of motion 

qk{t) = i^kQkiO) sinwfct +pfc(0) s'mujkt, (29) 

we find the correlation function 

{qk{t)qk{t')) = ujI cos UJkt cos UJkt'{qk{0)qk{0)) 
+ smu;ktsiiiuJkt{pk{0)pk{0)) 
= cos ujkit ~t'). (30) 

We may now assume that the oscillators qk{t) are the 
Fourier components a massless field, for instance the 
gravitational field whose frequencies are tOk — k, and 
whose random intial conditions are caused by the big 
bang. If the sum over k is simply a momentum integral, 
then ( |30| ) yields a white-noise correlation function (jj) for 

vit)- 

6. We have shown that it is possible to simulate the 



quantum-mechanical wave functions iplx^t) and the en- 
ergy spectrum of an arbitrary potential problem by clas- 
sical stochastic equations of motion, or by deterministic 
equations with random initial conditions. 

It remains to solve the open problem of finding a classi- 
cal origin of the second important ingredient of quantum 
theory: the theory of quantum measurement to be ex- 
tracted from the wave function tplx^t). Only then shall 
we understand how God throws his dice. 
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